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We use a nonlinear Schroedinger-Poisson equation to describe two interacting electrons with 
opposite spins confined in a parabolic potential, a quantum dot. We propose an effective form of 
the Poisson equation taking into account the dimensional mismatch of the two-dimensional electronic 
system and the three-dimensional electrostatics. The results agree with earlier numerical calculations 
performed in a large basis of two-body states and provide a simple model for continuous quantum- 
classical transition with increasing nonlinearity. Specific intriguing properties due to eigenstate 
non-orthogonality are emphasized. 
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I. INTRODUCTION 



Quantum-dots can be viewed as artificially structured 
atoms in heterojunctions or metal-oxide-semiconductor 
devices where few electrons are confined to a length 
comparable to the mesoscopic effective Bohr radius qb 
{as ^ 10~^ fim in the case of GaAs). Though the 
confinement can a priori occur in all three directions, 
some types of experimentally realized quantum-dots dis- 
play an extension in the x — y plane which is much 
larger than in the growth direction z of the underlying 
semiconductor structurei^i^ Therefore, these quantum- 
dots are usually regarded as artificial atoms with a disk- 
like shape. Since electron numbers N as low as one 
or two per dot have already been realized, "'^■'^ quantum- 
dot Helium consisting of two electrons trapped in the 
two-dimensional (2D) axisymmetrical harmonic potential 
V{r) — ^Muj^r^, where r"^ = + and M is the effec- 
tive electron mass, is actually considered as the simplest 
realistic model for an interacting quantum system^-i^ As 
itself or amongst other such few-electron systems, it has 
been extensively studied in the relationship with the de- 
velopment of nanotechnologiesi^'^'^ Both its exact 2D- 
3D analytical^ or 2D numerical^ solutions in the pres- 
ence of a perpendicular homogeneous external magnetic 
field are known, in particular by use of the separation 
of the Hamiltonian into its center-of-mass and relative- 
motion terms, due to the assumption of a parabolic par- 
ticle confinement. Oscillations between spin-singlet and 
spin-triplet ground states as a function of the magnetic 
field strength have been predicted2,i2, and experimentally 
observed?^ They are due to the interplay between the dot 
size and the strength of the magnetic field. Another im- 
portant competition occurs between the kinetic-energy 
matrix elements and the electron-electron Coulomb in- 
teraction ones when changing the characteristic length 
L of the quantum dot without changing its shape. In- 
deed, for small L, the Coulomb interaction becomes neg- 
ligible and the electrons behave like independent, un- 
correlated particles.^*' This happens in particular in the 
case of strong parabolic confinement u oo since then 
L^Iq = y/h/Muj 0. 



In this rich theoretical and experimental context, we 
wish to emphasize new physical results by use of a quite 
original - with respect to the above state of the art - 
differential approach based on the Schroedinger-Poisson 
(SP) definition of single-particle nonlinear eigenstates in 
quantum-dot Helium. The problem with such eigen- 
states is that, being the (stationary) solutions of the SP 
nonlinear differential system, they are not orthogonal. 
The whole matrix machinery of quantum mechanics then 
fails and we are left to return to its ab-initio fundamen- 
tal principles. In particular, the square scalar product 
V = (^al^b)^ ^ of two such nonlinear eigenstates 
and '^b defines the probability to find the system 
in either state when it is known to be in the other one 
(this probability is of course zero for orthogonal Hilber- 
tian eigenstates). Equivalently - and this will be pre- 
cisely shown below by use of the Fermi golden rule - 
, V yields the transition probability either from '^a to 
^'f, or reverse. Therefore, if '^a is, say, the fundamental 
eigenstate and \l/f, is an excited one, V either measures 
the probability of an absorption process (^'a ^'b) or 
an emission one (^j, '^a)- The following couple of 
questions are then addressed and tentative answers pro- 
vided: i) What energy is actually absorbed or emitted as 
a consequence of the non-orthogonality of the nonlinear 
eigenstates ^'a.fc? ii) Are the affiliated energy exchanges 
quantized and how? 

The present paper is built as follows. We first dis- 
play, by use of standard numerical tools, the remark- 
able properties of the solutions {u_^, CV} of the SP di- 
mensionless nonlinear differential system, where Af is the 
nonlinear control parameter whose value is given by the 
harmonic trap parabolicity uj. Then we numerically in- 
vestigate the square scalar product Vis = (ui l^a)^ ^ of 
its first two (zero-angular-momentum, for the sake of sim- 
plicity) eigenstates ui_3 when Af is increased. Subscripts 
always refer, in the present work, to those single-particle 
energy eigenvalues (in units of hcu) which correspond to 
the Af ^ linear limit. We find an interference- like 
pattern oc sin^(iA/'). Then we validate the above re- 
sults both by stressing the link between Vis and Fermi's 
golden rule for small nonlinearity Af < 1, as well as by 
displaying the transition of the system to the asymptotic 
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semi-classical Thomas- Fermi regime for high values of the 
nonlinearity A/" ^ 1. Finally, we test the reliability of the 
present SP description of quantum-dot Helium by com- 
parison with the existing numerical (Refs (jJl^ & 
respectively Figs [U [3l and [5]) and analytical (Ref. (@): 
Fig. lb) data concerning either the fundamental energy 
level or, like in Ref. its (electro)chemical potential. 
In all cases our nonlinear SP eigensolutions do agree sur- 
prisingly well (i.e. within the percent) with the existing 
exact results (for comparison, the Hartree and Hartree- 
Fock departures from the exact PGM fundamental en- 
ergy eigenvalue^ are respectively 44 % and 8 % in PGM's 
Fig. 3). 

Therefore we recover an important property that has 
already been emphasized in the N = 2 Coulomb case, 
both for atomic Helium^^ and hydrogen ion H~ M- 
Namely that the SP nonlinear differential description 
yields surprisingly accurate values for the ground state 
energy when compared to their corresponding mean-field 
Hartree-Fock ones. The reason seems to be the very par- 
ticular physical system that is actually constituted by a 
mere Cooper-like pair of opposite-spin electrons trapped 
in the same orbital bound state: one electron, say elec- 
tron a with orbital wave function ^'a, "feels" the repul- 
sive electrostatic potential that is being created by 
its fellow electron b with orbital wave function ^f^. This 
potential is defined by the classical Poisson equation 
V • V$h cx — I'I'bp while '^a is solution of the single- 
particle Schroedinger equation including both classical 
potentials, namely the external confining potential V{r) 
and the electrostatic potential The system is closed 
by the "bosonic orbital assumption" = \l/f,. As a con- 
sequence, there is no (positive) electron self-interaction 
energy contribution like in Hartree's mean-field descrip- 
tion. Neither does (negative) Hartree-Fock's exchange 
energy play a significant role, due to our S = opposite- 
spin assumption. Therefore, the only remaining dif- 
ference with respect to the exact corresponding energy 
eigenvalues might be due - or at least partially - to the 
next-order (negative) correlation effects. 



II. THE 2D RADIAL 
SCHROEDINGER-POISSON NONLINEAR 
ORDINARY DIFFERENTIAL SYSTEM 

The SP differential system results from coupling the 
single-particle stationary Schroedinger equation that de- 
fines the 2D orbital wave function ^'(a;, y) in the potential 

Vix,y) + ^{x,y) 



' 2M 



+ V{x, y) + $(a;, y) *(x, y) = /i'J'(x, y), (1) 



with the Poisson equation which solely defines the mu- 
tual electrostatic repulsive interaction ^(x,y) between 
the two particles^ii^ 



Since j^'p oc [length]"^, we must indeed introduce, for 
dimensional reasons related to Eq. ([2|), a characteristic 
energy which we wish to write as ^N'hcj. The corre- 
sponding dimensionless parameter Af = M{uj) will be 
defined below and is a typical measure of the SP non- 
linearity. It is important to keep in mind that the above 
system Eqs (Il][2]) is only relevant for particles in the same 
orbital state 5*. 

Assuming the 2D axisymmetrical parabolic confining 
potential V{x,y) = V{r) 



iMw^r^, we have: 



*(x, y) = *(r, (j>) = ^]{r) e 



imcj) 



(3) 



The wavefunction ^E* is thus the eigenstate of the angular- 
momentum operator —ihd/dcf) related to its eigenvalue 
mh. Its radial part ^^ir), which describes the 2D confine- 
ment of the electron system in the z = transverse plane 
with radial symmetry in agreement with the experimen- 
tal results, is defined by 
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(4) 



if we i ntroduc e the following dimensionless quantities 
(Jo = y^h/Muj is the characteristic parabolic length) 



X^V2-- u = 
lo 



in order to scale the Poisson nonlinearity (namely, the 
r.h.s. of Eq. ([2]) ) to unity, as evidenced by the r.h.s. of Eq. 
(|3)d). The dot stands for derivation with respect to the 
(dimensionless) radius X and, as already emphasized, the 
tilde superscript labels energy in units of fiuj. The single- 
particle probability of presence j^'p must be normalized 
to unity. Therefore M is the corresponding dimensionless 
norm of the solution u 
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XdX = TV. (6) 



III. 



THE NONLINEAR QUANTUM-CLASSICAL 
TRANSITION 



The dimensionless solution of our differential problem 
Eqs (HI yields u{r) and C{r) as functionals of M given by 
Eq. ([6]). It is defined by the initial conditions uq ~ u(0), 
iiQ = u{0), Co = C(0) and Co = (7(0). Amongst them, 
uo and Co are left free and will be chosen by numerical 
dichotomy, in order to yield regular bound-state eigen- 
states, i.e. u solutions defined by u{X) ^ for X ^ oo 
(for practical purposes, it will be sufficient to impose 
u{X) < IQ-'^ at X = 9: see Fig. H below). On the 
other hand, Co = (no potential cusp a,t X = 0) and the 
determination of the last remaining parameters, namely 
Mo, proceeds from step- by-step increase of both uq and 
Mo, starting from their linear limit where uq Vf^ ^ 1 



V2$(x,y) = -2nUhuj\^{x,y)f' 



(2) 



Mo < 1; uo 



■lin 



(0). 



(7) 
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Here 

u""(X)oce-*X^P„(X), (8) 

defines^ the Af ^ linear solutions in terms of the 
quantum numbers n and m and of the Laguerre poly- 
nomial Pn{X), namely Pq = 1; = 1 - ix^; = 
2 - + 1 jj^4. Pg ^ 6 - + - ... These 
single-particle linear parabolic states correspond to the 
energy eigenvalues 

E^"' = En,rn = (2n + |m| + 1) Uuj. (9) 

Figure [T] displays the Co versus uq spiraling trajectories 
for the two first m — eigenstates ui (circles) and 
(stars). Recall that the subscripts always refer in the 
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FIG. 1: (Color online). The convergence of the discrete 
SP nonlinear system towards the quasi-classical continuum 
Thomas-Fermi regime defined by the fixed point u = 1 and 
Co = lifX'^Jx^o = when Af increases from ~ 10~^ to ~ 10^ 
in the {Co vs uo} boundary condition phase space for the 
two first m = nonlinear eigenmodes defined by their cor- 
responding linear quantum numbers: namely (cf. Eq. (|9}) 
n = 0, m = (ui: circles) and n = 1, m = (113: stars). 

present work to the number of Hu quanta present in the 
Af ^ linear limit of the single-particle energy, in ac- 
cordance with Eq. [51 These trajectories in the {uq, Cq} 
plane are parametrized with respect to increasing values 
of the nonlinear parameter Af, i.e. with decreasing val- 
ues of the trap harmonicity uj. Indeed electron-electron 
interaction becomes relatively (with respect to quantum 
kinetic energy) more and more important when the two 
electrons are less and less confined (see above Part I). 
Actually Af varies in Fig. [T] from 10"^ (ug Q.l) to 
10^ [uq ^ 1) where one then reaches the quasi-classical 
asymptotic Thomas-Fermi regime. This regime is de- 
fined by neglecting the quantum kinetic derivative terms 
in Eq. (g^), thus yielding C{X) ~ X'^/i and hence 
Co = C(0) = 0, while uiX) EE 1 through Eq. (gla). There- 
fore the initial conditions for the two discrete modes ui 



and W3 converge towards the Thomas-Fermi fixed point 
{uq = 1; Co = 0} for A/" ^ 00 as evidenced by Fig. [1] 
Physically, this means that there is a continuous transi- 
tion, through the increase of nonlinearity in the system, 
from the Af < tt "pure" quantum regime where the quan- 
tum kinetic energy defined by the derivative terms in Eq. 
(H^) plays a major role towards the Af ^ tt classical one 
where the dimensionless Schroedinger equation Eq. (g^) 
reduces to its last-bracket classical-energy term. As a 
consequence, the A\f ^ 00 highly nonlinear case leads to 
the progressive merging of the two discrete energy lev- 
els ui and M3 into the single one whose initial conditions 
are defined by the fixed point displayed in Fig. [T] There- 
fore quantum eigenstate discreteness disappears, which is 
the hallmark of the classical regime: a continuous energy 
spectrum sets on about the uniform wavefunction profile 
u{X) = 1 and C{X) = 0, where the chemical potential 
equals the - here vanishing, due to our m = assumption 
- centrifugal potential plus the electrostatic interaction 
potential, as shown by Eq. ((S}:). 

The onset of the first corresponding oscillation in the am- 
plitude of the respective modes Ui{X) (continuous line) 
and us{X) (dotted line) is displayed in Fig. [21 
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FIG. 2: (Color online). Several nonlinear eigenstate profiles 
ui{X) (continuous line) and U3{X) (dashed line) for increas- 
ing values of the dimensionless quantum-dot size k, namely 
10~^ < k < 8.7, where k — lo/a* is the ratio of the char- 
acteristic harmonic length lo — y^h/Muj over the effective 
Bohr radius a* = /Me^ . The maximum-amplitude thresh- 
olds at Mo ~ 1.2 (resp. uo ~ 1.8) for the ground-state mode 
ui (resp. the excited mode us) and displayed by Fig. [1] are 
clearly visible (profile accumulation effect). 



IV. CHEMICAL POTENTIAL AND ENERGY: 
THE EXPLICIT DEFINITION Af{Lu) 

In the following, we shall only consider zero-angular- 
momentum m = states for the sake of simplicity (we 
have indeed checked that m 7^ nonlinear eigenstates are 
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equally well described by the above differential system: 
see below Fig. 2] where the quantum-dot spectra are dis- 
played versus their corresponding J\f for n < 3; m < 3). 
The SP virial energy E per particle corresponding to 
the nonlinear eigenstate u{X) is twice the expectation 
value (iMw^r^) of the external parabolic potential en- 
ergy (virial theorem for a harmonic potentiel). In terms 
of the dimensionless quantities defined in Eqs ([5]), it reads 



E 



1 



u^X^dX. 



(10) 



On the other hand, the chemical potential /i defined by 
Eq. ^ is that energy which is required in order to add 
the second electron to the single-electron quantum-dot 
(Koopman's theorem^). It can truly be regarded as the 
nonlinear eigenvalue of the SP differential system related 
to the corresponding nonlinear eigenstate ^' (or u in the 
reduced units defined in Eqs (O). Therefore we have 



/i = 2E-E 



lin 



(11) 



where i?''" — £"„.,„ is defined by Eq. [S] In the present 
work where we only consider m = eigenstates, the two 
first levels are Ei (resp. £^3) corresponding to n = (resp. 
n — 1) in units of tuv. The nonlinear integrodifferential 
system Eqs (|H ITT|) is closed by the use of Eq. ([St) at 
X = 0. This yields the (reduced) chemical potential 
/2 = ii/fiLU for m ~ 



^ = C(0) + $(0) = Co + — 



= Cn 



m 



dX, (12) 



where k = Iq/q* is the usual dimensionless dot size cor- 
responding to the harmonic length ~ y/h/Mu and 
the effective Bohr radius a* — /Me^ (ranging from 
a* = 67 nra for InSb to a* = 9.8 nm for GaAs). 

Equations (j4lfT2| self-consistently define, for any given 
value of the trap characteristic harmonic frequency lo (or, 
equivalently, its reduced size k), the solution u = Uaj{X), 
its norm Af = JV{uj) as well as its corresponding single- 
particle energy E = E{lo), together with the chemical 
potential jl = fl{uj). We numerically obtain, for instance 
for the ground state in the "quantum-regime" interval of 
values J\f < TT (see below) 



J\f{k) 



0.8839 fc 



0.4218 + 0.1247 /c 



(13) 



while its energy is 

E{JV) - 1-t- 0.24670 7V+ 0.03683 7V2-0.002177V^ (14) 

Therefore M{k) is a monotonic increasing function of the 
dot size, starting like N ^ 2k for small values of the dot 
size fc, while 



E E 



lin 



1 



M < 1 



k < 



1 



(15) 



In Ref. 'i, for instance, where Hlv = 3.37 meV for a 
GaAs parabolic quantum dot (M equals 0.067 electron 
mass while the charge is l/-\/12.4 electron charge), we 
have a* — 9.79 nm and Iq = 18.5 nm. Hence k = 
1.89. Then Eqs l[^3^^ respectively yield Af = 2.53 

and E — -Epcr particle = 1-83. Therefore iS'quantum dot — 

2i?pcr particle — 2(1.83) Huj — 12.33 mcV to be compared 
with PGM's value 12.28 meV: see Fig.[3]where the virial 
energy per particle (solid line), defined by Eq. (fTO|) . is 
plotted together with the Koopman one (dashed line), 
defined from Eqs ([9l fT2|) . The energy per particle E is 
defined by the intersection of both plots. Similarly, in 
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FIG. 3: (Color online). For GaAs quantum-dot Parahelium 
defined by the confinement hu = 3.37 meV, the dielectric 
constant of the bulk material k = 12.4 and the effective mass 
M = 0.067 me (where me is the electron mass),— we have 
A/" = 2.53 and E = 1.831 by use of (respectively) Eqs ^ 
& (|14|) . This last SP ground-state energy per particle value 
appears here as the intersection of its virial (continuous line) 
and its Koopman (dashed-dotted line) definitions as respec- 
tively provided by Eqs (fTOll & (|llj|12|l . PGM's exact numeri- 
cal value E = ^(12.28) meV /3.37 meV = 1.822 given in Ref. 
(5) is plotted as the horizontal segment. 

Ref. % hu ^ 2 meV yields k = 2.43. Hence TV = 2.97 
and E — 2.00, which yields fi = 2E — 1 = 3.00 in accor- 
dance with Eq. (fTT|) . and therefore /i = 3.00 ?ia; = 6.00 
meV which is in complete agreement with the Coulomb- 
interaction case (di — d2 — oo) of that reference. 

A remarkable property of the SP nonlinear eigenso- 
lutions is their "universal" limit behavior defined by 
Eq. for small A/", whatever the actual state's quan- 
tum numbers n and m in Eq. ([9]): see Fig. |4] where 
E = -I- \N is plotted in dashed line. Therefore the 
energy 
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(16) 



which was introduced for dimensional reasons into the 
Poisson equation ^ is simply the smallest additional 
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FIG. 4: (Color online). Illustration of Eq. CH) (dotted line) 
for the 1 < -E''" < 5 nonlinear eigenmodes (as defined by their 
linear A/" — > quantum numbers n & m: see Eq. (IS}). The 
stars display the greatest ito value reached for each mode: see 
Fig. 1 for the two first m — modes where the stars would 
respectively correspond to uo ~ 1.2 and ito ~ 1.7. From 
bottom to top (a & b superscripts label degeneracies which 
are lifted by nonlinearity) : {n — 0, m — 0}, {n = 0, m = 1}, 
{n = 1, m = 0}", {n = 0, m = 2}", {n = 1, m = 1}, {n = 



2, m = 0}^ {n = l,m = 2}" 



quantum-dot energy due to particle-particle interaction 
nonlinearity. Indeed Eq. (|15p yields for the 2E quantum- 
dot energy 



lim 2E 



2E 



lin 



(17) 



where i?'™ is the — linear energy per particle de- 
fined by Eq. Therefore A is the smallest interaction 
(or nonlinear) energy value in our two-electron SP sys- 
tem that comes in addition to the already existing "lin- 
ear" quanta hui, when Af ^ 0. To see whether A is a 
true "nonlinear quantum" of energy - i.e. whether the 
energy exchanges between the two levels Ei and E3 can 
be described in the terms of both Tico and A - demands 
to define the actual transition probability between these 
levels from the non-orthogonality of the corresponding 
eigenstates. This will be done in the next part. 



THE SCALAR PRODUCT {mlua) AND THE 
CORRESPONDING TRANSITION 
PROBABILITY 



together with 



/"OO /"OO 

/ ujX^ dX - Afi{l + Cf^) / ufdX 
Jo Jo 



ujdX 



(19) 



of both zero-angular-momentum eigenstates ui(X) and 
us^X) corresponding to those energies Ei and E3 which 
are respectively defined at A/1,3 ^ by £^0,0 and i?i,o in 
Eq. Equation (|19p is the matching condition 



V^J^ujdX V2j;^uldX 



(20) 



defined by Eqs ([Tl fT^ which states that the trap har- 
monicity w, or equivalently its quantum-dot dimension- 
less length k = y/h/Muj/lh^ /Me^), must be identical 
for the two modes ui,3 that enter the calculation of the 
scalar product (|18p . Practically, in the numerical simula- 
tions of Eqs (|18ll20p , we will consider Eq. (fT9| as verified 
if it is fulfilled within an 10~^ error. 

In order for the scalar product Eqs (|18ll20p to make 
physical sense, we wish to link it with standard time- 
independent linear perturbation theory in the case of 
small nonlinearity A/" <C 1, i.e. for small "perturba- 
tive" particle-particle interaction $ defined by Eqs ([T][2). 
Therefore we deduct that 



lim (wiIms)^ = Pi 3, 



where 



,lin 



^pcrt I 



4'") I 



sin^(iLJi3i) 



(21) 



(22) 



■'1 3 



together with huji^ = E3 — Ei ^ 2fku, yields Fermi 
golden rule's transition probability per particle Indeed 
both the energies Ei ^ -Ei(i_i) P^r particle {i — 1,3: 
cf. Eq. ([9])) and the corresponding normalized eigen- 
states Ui ~ m''" in Eq. (|22p are those corresponding to 
the unperturbed linear system (hence the superscript), 



namely^ m']'" = e ^ and 



,lin 



(1 



2£1 
2 



(cf. 



Eqs ([8])). The perturbation potential i/P°''* is equal 
to i$ (per particle: hence the factor i) where $ — 
i -f |>(3)j is the interaction potential that has been 
averaged over its two components and There- 
fore 



lin\ 



fVport 
^13 



^[^13] 



(23) 



Let us define the (normalized) scalar product 
1 



\/A/iA/3 Jo 



U1U3 XdX, 



(18) 



where the matrix elements $1*3 (i = 1, 3) have been cal- 
culated by use of the above-mentioned linear normalized 
eigenstates m'i™. In the stationary perturbative regime 
where W13 i ~ 2u!t ^ 1 3> HP"^t/h, the time-dependent 
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term in Eq. ([22)l can be replaced by its averaged value ^ , 

yielding T'ls = M^ff ]' = M^ii + ^ilT- Therefore 
Eq. (pTj) becomes 

^iin^li3 = I [$W + 1.^3)] = 2V2 1(^.1 1^3)1 . (24) 

Figure El displays the r.h.s. of Eq. ([24]) (continuous line) 
versus its l.h.s. (dashed-dotted line) and shows the nu- 
merical fulfi-Uement of this condition. 
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FIG. 5: (Color online), the convergence for small nonlinearity 
Af ~ A/i A/a < 1 of the SP nonlinear model towards Fermi's 
"golden rule" probability transition defined by Eqs (|18I22|I , as 
displayed by Eq. (|24|| . 

Therefore the square scalar product defined by Eqs 
(|18ll20p yields, in the limit of small nonlinearity A/1,3, the 
transition probability Vi 3 from the fundamental nonlin- 
ear eigenstate ui to the excited one U3 or reverse. On 
the other hand, we showed in Section III that the quasi- 
classical Thomas-Fermi regime yields \mij\f^oo "1,3 = 1 
(see Fig. [ij. The two modes M1.3 then become equiv- 
alent. Consequently the transition probability between 
them should obviously become equal to unity, which is 
consistent with [(ui |u3)^]ki^ii3 = 1 from definition 
Therefore it seems quite natural to extrapolate to all 
values of the nonlinearity A/1,3 the physical meaning of 
(mi|m3)^ in terms of the transition probability V13 as de- 
fined by Eqs (^111^ . 

VI. QUANTUM TRANSITIONS BETWEEN 
TWO m = NONLINEAR EIGENSTATES 

Let us now proceed to the investigation of the quantum 
transitions between the two nonlinear eigenstates ui,3 by 
use of the numerical calculation of the scalar product 
defined by Eqs (|18ll20p . It consists in increasing the non- 
linearity through a three-loop iterative scheme from the 
A/" <C 1 linear regime. The two first loops define each 
eigenstate ui,3 which vanish with a 10"'' accuracy at 



X 9 which is our numerical value for X 00 (see 
Fig. [5]) while the third one evaluates the matching condi- 
tion Eq. ([TO]) within 10~^ and then calculates the scalar 
product given by Eq. (|18p . The integrals which appear 
in Eqs (|18II20[) are transformed into additional first-order 
ordinary differential equations with vanishing initial con- 
ditions whose solutions are taken at X ~ 9. Then the 
whole resulting differential system is numerically inte- 
grated by use of standard tools-^i 

Figure[6]displays the following remarkable interference- 
like pattern with respect to the ground-state nonlinear 
parameter A/i. Intriguing enough, since the present SP 
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FIG. 6: (Color online). Equation (|25p 's square scalar product 
{Mi|M3)'^/7ra normalized by ttq — 7r/137.036, as compared 
with its sin^ approximation (dotted line), together with the 
nonlinear resonance condition Eq. (|26[) (in "tilde" units of 
fvjj): i) intersection of the two upper plots that yields /is — 
\i\ ~ Ai ~ 7i"/2, and ii) intersection of the dashed-dotted 
lower plot with zero that yields Ai ~ A3, i.e. M\ ^ Nz, at 
the resonance nonlinearity A/i ^ tt. 

differential model is non-relativistic (there is no velocity 
of light in it), it is best scaled by use of the numerical 
value of the fine-structure-constant a. = /he multiplied 
by TT, namely 7r/137.036 = 2.2925... IQ-^ 

— luAuzf - 1.0005 sin^rO.SOeOAAil. (25) 

TTQ! 

We note that, when Mi > 4, the departure from the r.h.s. 
of Eq. (|25|) (dotted line in Fig. [6]) becomes significant 
as the transition toward the asymptotic quasi-classical 
Thomas- Fermi regime sets on. On the other hand. Fig. [5] 
displays the following chemical-potential gap transition 
process 

/i3-/xi~A^, (26) 

TT 

where A^r ~ Ai ^ A3 ^ — ?ii^7r is the common charac- 
teristic energy A, defined by Eq. of the two eigen- 
states ui,3 about the maximum (of amplitude 1.0005 Tra) 
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of their square scalar product (wijus)^, i.e. at the very 
pecuHar quantum-dot nonhnearity Ni ^ As ~ tt re- 
lated to the particular oj — uJt^ trap parabolicity. This 
value corresponds to the specific parabolic confinement 
/iw^ ~ 0.14e where e — Me'^ /h^ is the effective quantum- 
dot's atomic energy unit: e = 11.86 meV for AsGa, thus 
yielding Tiw^r ^ 1-66 meV and A^r ^ 2.61 meV while 
e = 27.21 eV if the dielectric constant of the bulk material 
equals unity, then yielding TilUt^ ~ 3.80 eV and therefore 
- 5.97 eV. According to Eq. HIl), Eq. dH]) yields the 
corresponding quantization rule for the 2E quantum-dot 
energy at uj ^ lOt^ 

2iE3-E,)^2hLu^ + A^. (27) 

Equations (fT7| and ([27|) show that the characteristic en- 
ergy A which scales the electrostatic particle-particle in- 
teraction through the nonlinear differential Poisson equa- 
tion ([2]) is in fact a true "nonlinear quantum". In- 
deed, on the one hand, it is the smallest particle-particle 
interaction energy present in the system at vanishing 
nonhnearity A/" — > 0. On the other hand, the maxi- 
mum of the (miIms)^ transition probability between the 
two states is reached at resonance, i.e. either when A 
equals their nonlinear-eigenvalue chemical-potential gap 
or when their quantum-dot energy gap is but the mere 
sum of the two standard "linear" radial quanta htu and 
A. 



VII. CONCLUSION AND PERSPECTIVES 

In the present paper, we have described the parabolic 
quantum dot by use of the nonlinear differential eigen- 
problem Eqs ([TJS]) and emphasized its relevance with re- 
spect to all existing corresponding results in the liter- 
ature. This Schroedinger-Poisson (SP) differential sys- 
tem yields new quantum concepts such as the non- 
orthogonal nonlinear eigenstates ^ and their correspond- 
ing chemical-potential nonlinear eigenvalues /xJ^ In or- 
der to comply with the dimensional self-consistence be- 
tween the two-dimensional electronic system and its 
three-dimensional electrostatics, we scaled the Poisson 
equation according to the characteristic energy A — 
■^AfTiLU where A/" is a normalized (see Eq. measure 
of the system nonhnearity. We showed that A is actu- 



ally the true "nonlinear energy quantum" of the system 
for: i) it is the smallest additional "nonlinear" particle- 
particle interaction energy with respect to the standard 
"linear" radial harmonic quantum huj when Af —^ 
(see Eq. ([TT)) ): ii) it fits with that nonlinear-eigenvalue 
(or chemical-potential) gap between the two first zero- 
angular-momentum eigenstates which occurs about the 
maximum of their square scalar product (Mijws)'^ (see 
Eqs (|26II27|) ). i.e. about the maximum of their transition 
probability V (as a consequence of Fermi's golden rule). 

Further developments of the present work should (non 
exhaustively) address the two following experimental, nu- 
merical as well as theoretical topics: 

1) Could the nonlinear resonance defined by Eqs PBI ^ 
I27p at the very particular trap parabolicity htu — hajj^ 
{— 1.66 meV for GaAS) be observable and how? It would 
be a definite plus for the present model to provide an 
opportunity for experimental verification. 

2) The na scaling adopted in Eq. (pS)) seems extremely 
accurate. Indeed the square scalar product maximum 
divided by n approaches the numerical value 1/137.036 
of the fine-structure-constant /he within 0.05 % in the 
latest state of our numerical simulations 

^(ui|.3)L. = Yi7- 

It seems hard to believe that Eq. (|28p is but the re- 
sult of a mere numerical coincidence. Rather, we wish 
to point out that Eq. ([^5)1 might echo Feynman's em- 
phasis of such a "magic number" fii This stunning non- 
relativistic property of the nonlinear-eigenstate square- 
scalar-product scaling will be further investigated in a 
future publication. 
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